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Abstract 

Regge calculus is a powerful method to approximate a continuous manifold 

by a simplicial lattice, keeping the connectivities of the underlying lattice 

fixed and taking the edge lengths as degrees of freedom. The discrete Regge 

model employed in this work limits the choice of the link lengths to a finite 

number. To get more precise insight into the behavior of the four-dimensional 

discrete Regge model, we coupled spins to the fluctuating manifolds. We 

examined the phase transition of the spin system and the associated critical 

exponents. The results are obtained from finite-size scaling analyses of Monte 

Carlo simulations. We find consistency with the mean-field theory of the Ising 

model on a static four-dimensional lattice. 
PACS: 04.60.Nc, 05.50.+q 
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I. INTRODUCTION 



Spin systems coupled to fluctuating manifolds are analyzed as a simple example for mat- 
ter fields coupled to Euclidean quantum gravity. The gravitational action is unbounded 
from below due to conformal fluctuations. But that does not necessarily render its quan- 
tum theory useless or the path integral undefined. Indeed, the existence of a well defined 
phase was the first and probably most important result of the numerical simulations in four 
dimensions . Its existence and stability were explored in some detail using the standard 
Regge calculus with continuous link lengths. It turned out that the well defined phase is 
stable against variations of the measure and the lattice size |§. 

In the discrete Regge model the problem of an unbounded action is not present as in 
standard Regge calculus. Because of the restriction of possible quadratic link lengths to two 
values |4]|| in the discrete Regge model the action can only reach an extreme but finite value. 
The expectation values do not diverge if the well defined phase is left. What happens is 
that the lattice "freezes" at large positive and negative values of the gravitational coupling, 
as expected for a spin system ||. To get more precise ideas about the behavior of the four- 
dimensional discrete Regge model, we coupled Ising spins to the fluctuating manifolds. We 
examined the phase transition of the spin system and the associated critical exponents. 

The rest of the paper is organized as follows. In Sec. |T| we introduce the discrete Regge 
model and give some details of the analyzed observables. The results of the Monte Carlo 



simulations are presented in Sec. Oil and concluding remarks can be found in Sec. IV 



II. MODELS AND OBSERVABLES 

The situation for the discrete Regge model is both structurally and computationally much 
simpler than the standard Regge calculus with continuous link lengths. The restriction of 
the edge lengths to just two values was carefully examined in 2D where an interpolation from 
Z<i to was performed J7j. It turned out that the phase transition with respect to the 
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cosmological constant is universal. This was tested for pure gravity in 2D and is expected 
to be the case also in 4D. Compared with standard Regge calculus, numerical simulations 
of the discrete Regge model can be done more efficiently by implementing look-up tables 
and using the heat-bath algorithm. In the actual computations we took the squared link 
lengths as q^ = q\ = 6/(1 + eai) with er z ± 1. The Euclidean triangle inequalities are 
satisfied automatically as long as e < e max . Because a four-dimensional Regge skeleton with 
equilateral simplices cannot be embedded in flat space, bi takes different values depending on 
the type of the edge /. In particular bi = 1, 2, 3, 4 for edges, face diagonals, body diagonals, 
and the hyperbody diagonal of a hypercube. 
We investigated the partition function 

Z = E / D[q] exp [-l(q) - KE(q, s)] , (2.1) 
W 

where I(q) is the gravitational action, 

/(<?) = -A,E^ + A E^- ( 2 - 2 ) 

t i 

The first sum runs over all products of triangle areas A t times corresponding deficit angles 
5 t weighted by the gravitational coupling (3 g . The second sum extends over the volumes Vi 
of the 4-simplices of the lattice and allows one together with the cosmological constant A to 
set an overall scale in the action. The energy of Ising spins Si G Z2, 

E( q ,s) = W^ ( -^^, (2.3) 

z (ij) qi i 



is defined as in two dimensions [pHlO|l, with the barycentric area Ay associated with a link /j,-, 



Mi = E \ A f (2-4) 

t ZD lij 

We chose the simple uniform measure as in the pure gravity simulations [H]: 

0[g] = n^Oft)- (2-5) 
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The function T ensures that only Euclidean link configurations are taken into account, i.e., 
T — 1 if the Euclidean triangle inequalities are fulfilled and T = otherwise. This is always 
guaranteed for the discrete Regge model by construction. 

For every Monte Carlo simulation run we recorded the time series of the energy density 
e = E/N and the magnetization density m = J2i s i/No, with the lattice size N = L 4 . To 
obtain results for the various observables O at values of the spin coupling K in an interval 
around the simulation point Kq, we applied the reweighting method ||11|| . Since we recorded 
the time series, this amounts to computing 



\Ko 

with AK = K - K . 



(0)\k= tZZakexi > ( 2 - 6 ) 



With the help of the time series we can compute the specific heat, 

C(K)=K 2 N ((e 2 )-(e) 2 ) , (2.7) 

the (finite lattice) susceptibility, 

X (K) = iV ((m 2 ) - (H> 2 ) , (2.8) 

and various derivatives of the magnetization, d(\m\)/dK, dln(\m\) / dK , and dln{m 2 )/dK. 
All these quantities exhibit in the infinite-volume limit singularities at K c which are shifted 
and rounded in finite systems. We further analyzed the Binder parameter, 

Ul{K) = 1-\^%. (2.9) 

It is well known that the Ul(K) curves for different lattice sizes L cross around (K C ,U*). 
This allows an almost unbiased estimate of the critical spin coupling K c . 

III. SIMULATION RESULTS 



In four dimensions, after initial discussions [I^-IjJ it is generally accepted that the 
critical properties of the Ising model on a static lattice are given by mean-field theory with 
logarithmic corrections. The finite-size scaling (FSS) formulas can be written as |15|.[T( 



£cxL(logL)4, (3.1) 

X ocL 2 (logL)^ =(L(logL)3)^, (3.2) 

Coc(logL)i, (3.3) 

K c (oo) - K C (L) oc L- 2 (logL)-^ = (L(logL)^)- 1 ^ , (3.4) 

where the critical exponents of mean-field theory are a — 0, (3 — 1/2, 7 = 1, and v = 1/2. 
To get more precise ideas about these logarithmic corrections, we first simulated the four- 
dimensional Ising model on a regular lattice. After this comparative study we turned to the 
four-dimensional discrete Regge model with spin fields. 

A. Ising spins on a regular lattice 

We studied the four-dimensional Ising model on a hypercubic lattice with linear size 
L = 3 - 16, 18, 20, 24, 28, 32, 36, 40, using the single-cluster update algorithm (Wolff) |7|. 
The simulations were performed with the knowledge of the value of the critical temperature 



obtained in previous Monte Carlo simulations and high-temperature series analyses |18 



K c = —— = 0.149 694 ±0.000 002. (3.5) 
k B T c 

We performed n(L) oc N Q /(C(L)) cluster updates between measurements for lattices 
L < 24, with the averaged cluster size (C(L)). After an initial equilibration time of about 
1 000 xn(L) cluster updates we took about 50 000 measurements for each of the small lattices. 
For the larger lattices we measured after each cluster update, therefore, we took about 
500 000 measurements after an initial equilibration of 100 000 cluster updates. Analyzing the 
time series we found integrated autocorrelation times for the energy and the magnetization 
in the range of unity for the small lattices L < 24 and in the range of (4 — 8) x L for the 
larger lattice sizes. The statistical errors were obtained by the standard jack-knife method 
using 50 blocks. 

Applying the reweighting technique we first determined the maxima of C, x, d{\m\) / dK , 
dln(\m\)/dK, and dln(m 2 ) /dK. The locations of the maxima provide us with five se- 



quences of pseudo-transition points K mayL (L) for which the scaling variable x = (K c — 
K max {L)){L{\ogL)T2) 1 / u should be constant. Using this fact we then have several possi- 
bilities to extract the critical exponent v from (linear) least-square fits of the FSS ansatz 
with multiplicative logarithmic corrections considering Eq. 



dU L /dK={L{\ogL)-^Ux) or (3.6) 
d\n(\m\P)/dK - (L(logL)^f p {x), (3.7) 

to the data at the various K m3fX (L) sequences. For comparison we also performed fits of a 
naive power-law FSS ansatz 

dU L /dK = L 1/u f (x) or (3.8) 
dln(\m\ p )/dK L 1/v f p {x). (3.9) 

The exponents 1/u resulting from fits using the data for the N largest lattice sizes are 
collected in Tables I and II. Q denotes the standard goodness-of-fit parameter. For all 
exponent estimates the FSS ansatz with the logarithmic corrections leads to the weighted 
average 1/u — 1.993(3), which is in perfect agreement with the mean-field value 1/u = 2, see 
Fig. |T] (a). With the naive power-law ansatz one also gets an estimate for 1/u in the vicinity 
of the mean- field value, but this is clearly separated from the mean-field result, verifying 
the significance of the multiplicative logarithmic correction, cf. Table II. 

Assuming therefore u = 0.5 we can obtain estimates for K c from linear least-square fits 
to the scaling behavior of the various K max sequences, as shown in Fig. || (b). Using the 
fits with L > 12, the combined estimate from the five sequences leads to K c = 0.149 697(2), 
which is in agreement with the results using Monte Carlo computer simulations [18[ and 
series expansions [|1£],^(J. 



Knowing the critical coupling we may reconfirm our estimates of 1/u by evaluating the 
above quantities at K c . As can be seen in Tables I and II, the statistical errors of the FSS 
fits at K c are similar to those using the K max sequences. However, the uncertainty in the 
estimate of K c has also to be taken into account. This error is computed by repeating the 
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fits at K c ± AK C and indicated in Tables I and II by the numbers in square brackets. In the 
computation of the weighted average we assume the two types of errors to be independent. 
As a result of this combined analysis we obtain strong evidence that the exponent v agrees 
with the mean- field value of v = 0.5. 

To extract the critical exponent ratio 7/V we use the scaling 

Xma x = (£(l0gL)3)^ (3.10) 

as well as the scaling of x at K c , yielding in the range L = 10—40 estimates of 7/V = 2.037(9) 
with Q = 0.95 and 7/V — 2.008(5) [5] (Q = 0.46), respectively. These estimates for 7/V 
are consistent with the mean-field value of 7/V = 2. In Fig. |2] (a) this is demonstrated 
graphically by comparing the scaling of Xmax with a constrained one-parameter fit of the 
form Xmax = c(L(logL)3) 2 with c = 0.526(2) (Q = 0.38, L > 6). 

Concerning the specific heat we expect in the case of the mean-field exponent a = a 
logarithmic divergence of the form 

C(x,L) = A(x) + B(x)(\ogL)K (3.11) 

Indeed, the data at the different fixed values of x can all be fitted nicely with this ansatz. 
In particular, for the fit of C max with 17 points (L > 6) we obtain A = —0.324(32), B = 
1.038(23), with a total \ 2 — H-7(Q = 0.70). We also tried an unbiased three-parameter fit 
using the ansatz 

C{x,L) = A'{x)+B'{x){\ogL)^ x \ (3.12) 

which in the case of the fit of C max and 16 data points gives A' pa —0.36, B' ~ 1.75, and 
k = 0.33(40), with a slightly improved total y 2 = 10. 8 (Q = 0.62). In Fig. |] (b) we compare 
these two linear least-square fits. It should be noted, however, that the three-parameter fit is 
highly unstable and exhibits strong correlations between the three parameters. To illustrate 
this instability we plot in Fig. |3] (a) the exponent k as a function of the smallest lattice size 
Limn, being the lower bound of the fit range [L min , 40]. For comparison we show in Fig. [3] (b) 
the behavior of 7/1/ results of the fit corresponding to Eq. (|3.10| ). 
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fit type 


N 


l/i/ 


Q 


dU/dK at Kg^ 
d\n(\m\)/dK at xS m|> 
d\n(m 2 )/dK at 
weighted average 


22 
18 
18 


1.990(5) 
1.989(4) 
1.998(5) 
1.993(3) 


0.68 
0.94 
0.93 


dU/dK at K c 
d\n{\m\)/dK at K c 
d\n(m 2 )/dK at K c 
weighted average 


18 
18 
18 


1.991(5)[1] 
1.992(5) [2] 
2.000(5) [2] 
1.995(3) 


0.73 
0.94 
0.94 


overall average 




1.994(2) 





Table I. Fit results for \ jv with a power-law ansatz with logarithmic corrections, using 
the data for the N largest lattices. The average is computed by weighting each entry with its 
inverse squared error. For the fits at our best estimate of K c = 0.149 697(2) the uncertainty 
due to the error in K c is indicated by the numbers in square brackets. 



fit type 


N 


\jv 


Q 


dU/dK at Kg ax 
d\n(\m\)/dK at^;° (H> 
d\n(m 2 )/dK at 
weighted average 


18 
18 
18 


2.041(9) 
2.056(5) 
2.066(5) 
2.059(3) 


0.81 
0.75 
0.63 


dU/dK at K c 
d\n(\m\)/dK at K c 
d\n(m 2 )/dK at K c 
weighted average 


18 
18 
18 


2.043(9) [2] 
2.059(5) [2] 
2.067(5) [2] 
2.061(4) 


0.81 
0.81 
0.70 


overall average 




2.060(2) 





Table II. Fit results for l/v with a pure power-law ansatz using K c = 0.149 697(2). The 
averages and statistical errors are computed as in Table I. 
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FIG. 1. (a) Least-square fits of the FSS ansatz with logarithmic corrections at the maxima 
locations. Together with the fits at K c this leads to an overall critical exponent 1/u = 1.994(2). 
(b) FSS extrapolations of pseudo-transition points if max vs. (L(logL)T2)~ 1 / ly , assuming v = 0.5. 
The error- weighted average of extrapolations to infinite size yields K c = 0.149 697(2). 



10000 



1000 



100 



10 



a) 



y/v = 2 



10 
L 



CO 

E 
O 



1.5 
1.4 
1.3 
1.2 
1.1 
1 

0.9 
0.8 
0.7 
0.6 
0.5 



J? 



b) 



A + B (log L) 
A + B (log L) 1 



1/3 



10 15 20 25 30 35 40 45 
L 



FIG. 2. (a) FSS of the susceptibility maxima Xmax- The exponent entering the curve is set 
to the mean-field value j/u = 2 for regular static lattices, (b) FSS of the specific-heat maxima 
C max - The logarithmic fit C max = A + B(logL) K and a constrained logarithmic fit assuming the 
mean-field prediction n = 1/3 are almost indistinguishable on the scale of the figure. 
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FIG. 3. (a) Instability of the logarithmic fit C max = A + B(\ogL) K . The exponent k as a 
function of the lower bound of the fit range is plotted, (b) y/v as a function of the lower bound of 
the fit range is plotted for the fit Xmax oc {L(\og L) i~yi/ v . 



B. Ising spins on a discrete Regge model 

We simulated the gravitational degrees of freedom (the squared link lengths) of the 
partition function ( |2.1|) using the heat-bath algorithm. For the Ising spins we employed 
again the single-cluster algorithm. Between measurements we performed n — 10 Monte 
Carlo steps consisting of one lattice sweep to update the squared link lengths followed 
by two single-cluster flips to update the spins Sj. 

The simulations were done for e = 0.0875, cosmological constant A = and two different 
gravitational couplings, (3 9 = —4.665 and j3 9 = 22.3. These two /3 9 -values correspond to 
the two phase transitions of the pure discrete Regge model ||, as shown in Fig. |J The 
transition at positive gravitational coupling of the standard Regge calculus was described 



in great detail in Ref. |2]] and shown to be of first order whereas the transition at negative 
coupling turned out to be of second order for the discrete Regge model || . Together with an 
eventual second-order transition of the Ising part, the latter one is a candidate for a possible 
continuum limit. The lattice topology is given by triangulated tori of size iVo = L 4 with 
L = 3 up to 10. From short test runs we estimated the location of the phase transition of 
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the spin model and set the spin coupling Kq = 0.024 rs K c in the long runs for both values 
of P g , see Fig. § 

After an initial equilibration time we took about 100 000 measurements for each lattice 
size. Analyzing the time series we found integrated autocorrelation times for the energy and 
the magnetization in the range of unity for all considered lattice sizes. As in the simulations 
of the regular lattices the statistical errors were obtained by the standard jack-knife method 
using 50 blocks. 

Completely analogously to the Ising system on a regular lattice we applied reweighting 
to locate the maxima and used the FSS formulas (3.6)- (|3.9|) . The exponents 1/v resulting 



from fits using the data for the N largest lattice sizes are collected in Tables III and IV for 
(3 g = —4.665, and in Tables V and VI for f3 g = 22.3, respectively. For the simulations at 
(3 g = —4.665 all exponent estimates with the logarithmic corrections and consequently also 
their weighted average 1/v = 2.025(6) are in agreement with the mean-field value 1/v = 2, 
see Fig. |6| (a). For (3 g = 22.3 the scatter in the estimates is similar and the weighted average 
1/v = 2.028(6) is again compatible with 1/v = 2. With the naive power-law ansatz one also 
gets an estimate for 1/v in the vicinity of the mean-field value, but this is clearly separated 
from the mean- field result, cf. Tables IV and VI. 

Assuming therefore v = 0.5 we can obtain estimates for K c from linear least-square 
fits to the scaling behavior of the various if max sequences, as shown in hboxFig. || (b) for 
(3 g = —4.665. Using the fits with L > 4, the combined estimate from the five sequences 
leads to K c = 0.02464(4) for (3 g = -4.665, and for (3 g = 22.3 we find K c = 0.02339(4), again 
with L > 4. 

With the knowledge of the critical couplings we may reconfirm our estimates of 1/v by 
evaluating the above quantities at K c . As can be inspected in Tables III and V, we obtain 
from this combined analysis strong evidence that the exponent v agrees with the mean-field 
value of v — 0.5. 

To extract the critical exponent ratio j/v we use the FSS formula fl3.10p for Xmax as 



well as the scaling of x & t K c , yielding for (3 g = —4.665 in the range L = 4 — 10 estimates 

11 



of 7/z/ = 2.039(9) with Q = 0.42 and 7/1/ = 2.036(7) [4] (Q = 0.85), respectively. The 
corresponding values for (3 g = 22.3, using the same fit range, are 7/1/ = 2.052(8) (Q = 0.57) 
and 7/z/ = 2.052(6) [4] (Q = 0.01). These estimates for 7/z/ are compatible with the mean- 
field value of 7/z/ = 2. In Fig. [7| (a) this is demonstrated by comparing the scaling of Xmax 
with a constrained one-parameter fit of the form Xmax — c(L(log L)^) 2 with c = 4.006(10) 
(Q = 0.17, L > 6) for /3 g = -4.665 and c = 4.244(10) (Q = 0.001, L > 6) for /3 g = 22.3, 
respectively. 

The data for the specific heat C at the critical spin coupling K c are presented in 
Fig. [7| (b). The fact that C increases very slowly with the size of the lattice means that 
one will need data from bigger lattices and more statistical accuracy to get an estimate or 
bound for the critical exponent a from a direct measurement of C . Especially, if we assume 
a logarithmic divergence of C as in the four-dimensional Ising model on regular lattices, we 
need lattices of comparable size, cf. Fig. |] (b). 
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FIG. 4. Expectation values of the average squared link lengths as a function of the gravitational 
coupling for the pure discrete Regge model on a 4 4 lattice. 
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FIG. 6. (a) Least-square fits of the FSS ansatz with logarithmic corrections for (3 g = —4.665 
lead to an overall critical exponent \jv = 2.025(4). (b) FSS extrapolations of pseudo-transition 
points K max vs. (L(log L)T2 )~ l / u for j3 9 = —4.665, assuming u = 0.5. The error-weighted average 
of extrapolations to infinite size yields K c = 0.02464(4). 
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fit type 


N 


l/i/ 


Q 


dU/dK at Kg^ 
<M{\m\)/dK at K$ H) 
d\n(m 2 )/dK at 
weighted average 


8 
7 
7 


2.003(10) 
2.032(10) 
2.038(10) 
2.025(6) 


0.47 
0.59 
0.55 


dU/dK at K c 
<M(\m\) / dK at K c 
d\n{m 2 )/dK at K c 
weighted average 


7 
7 
7 


1.981(17)[13] 
2.027(9) [2] 
2.034(9) [2] 
2.025(6) 


0.70 
0.95 
0.85 


overall average 




2.025(4) 





Table III. Fit results for \jv with a power-law ansatz with logarithmic corrections for 
P g = —4.665, using the data for the N largest lattices. The average is computed by weighting 
each entry with its inverse squared error. For the fits at our best estimate of K c = 0.02464(4) 
the uncertainty due to the error in K c is indicated by the numbers in square brackets. 



fit type 


N 


l/i/ 


Q 


dU/dK at iCx 
d\n{\m\)/dK at K^ ml) 
d\n(m 2 )/dK at 
weighted average 


7 
7 
7 


2.068(18) 
2.122(10) 
2.128(10) 
2.118(7) 


0.60 
0.37 
0.35 


dU/dK at K c 
<M(\m\)/dK at K c 
d\n(m 2 )/dK at K c 
weighted average 


7 
7 
7 


2.068(18)[12] 
2.116(9)[2] 
2. 124(9) [2] 
2.116(7) 


0.59 
0.83 
0.64 


overall average 




2.117(5) 





Table IV. Fit results for l/v with a pure power-law ansatz for f3 g = —4.665. The averages 
and statistical errors are computed as in Table III. 
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fit type 


N 


l/i/ 


Q 


dU/dK at Kg^ 
d\n(\m\)/dK at K$ H) 
d\n(m 2 )/dK at 
weighted average 


8 
7 
7 


1.981(10) 
2.043(9) 
2.049(9) 
2.028(6) 


0.64 
0.61 
0.67 


dU/dK at K c 
<M(\m\) / dK at K c 
d\n(m 2 )/dK at K c 
weighted average 


8 
7 
7 


1.993(10)[1] 
2.039(9) [2] 
2.045(9) [2] 
2.027(6) 


0.76 
0.32 
0.49 


overall average 




2.028(4) 





Table V. Fit results for \jv with a power-law ansatz with logarithmic corrections for 
P g = 22.3. The average is computed by weighting each entry with its inverse squared error. 
For the fits at our best estimate of K c = 0.02339(4) the uncertainty due to the error in K c 
is indicated by the numbers in square brackets. 



fit type 


N 


l/i/ 


Q 


dU/dK at Kg^ 
d\n(\m\)/dK at K^ ml) 
d\n(m 2 )/dK at 
weighted average 


7 
7 
7 


2.086(11) 
2.134(10) 
2.139(9) 
2.122(6) 


0.72 
0.54 
0.59 


dU/dK at K c 
<Sn{\m\)/dK at K c 
d\n(m 2 }/dK at K c 
weighted average 


8 
7 
7 


2.098(10)[1] 
2.130(9)[2] 
2.135(9)[2] 
2.122(6) 


0.57 
0.35 
0.48 


overall average 




2.122(4) 





Table VI. Fit results for \ jv with a pure power-law ansatz for /3 g = 22.3. The averages 
and statistical errors are computed as in Table V. 
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FIG. 7. (a) FSS of the susceptibility maxima Xmax- The exponent entering the curve is set 
to the mean-field value 7/1/ = 2 for regular static lattices, (b) Specific heat at the critical spin 
coupling as a function of the lattice size L. 
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IV. CONCLUSIONS 



We have performed a study of the Ising model coupled to fluctuating manifolds via Regge 
calculus. Analyzing the discrete Regge model with two permissible edge lengths it turns out 
that the Ising transition shows the predicted logarithmic corrections to the mean-field theory. 
The critical exponents of the phase transition of the Ising spins on a static lattice as well as 
on a discrete Regge skeleton P2] are consistent with the exponents of the mean-field theory, 
a = 0, /3=|,7 = 1, and v — ~. In summary, our consistent analysis with uniform computer 
codes yields that the phase transition of a spin system coupled to a discrete Regge skeleton 
exhibits the same critical exponents and the same logarithmic corrections [l(J as on a static 
lattice. 
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